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We report a study of B — ¥ (J/ipj)K and B —¥ (ip'^jK decay modes using 772 x 10 6 BB events 
collected at the T(4S) resonance with the Belle detector at the KEKB energy-asymmetric e + e~ 
collider. We observe X(3872) — > J/ip'y and report the first evidence for \c2 — > J/tpy in B — > (Xcc'y)K 
decays, while in a search for X(3872) — > tp'"/ no significant signal is found. We measure the branching 
fractions, B(B ± X(3872)if i )S(X(3872) -> J/i/rf) = (f .78tg'jt±0.12)xKT 6 , B{B ± -> x^K ± )= 
(1.11±°J| ± 0.09) x f0" 5 , B(B ± -> X(3S72)K ± )B(X (3872) ->• ^'7) < 3.45 x 10~ 6 (upper limit at 
90% C.L.) and also provide upper limits for other searches. 

PACS numbers: 13.20.Gd, 13. 20. He, 14.40. Gx 



The A (3872) state was observed by the Belle Collabo- 
ration Q in 2003, and later confirmed by CDF @, D0 @ 
and BaBar The fact that it was not seen in decays 
to Xci7i Xc27, and J/iprj final states suggests that the 
A (3872) is not a conventional qq meson state that can be 
explained by a simple quark model Because of its 

narrow width and the proximity of its mass, 3871.5 ± 0.2 
MeV/c 2 to the D*°TP threshold, the A(3872) is a 
good candidate for a DD* molecule Q. Other possi- 
bilities have also been proposed for the A (3872) state, 
such as tetraquark ccg hybrid meson [l(| and vector 
glueball models (llj . 

Radiative decays of the A (3872) are important in un- 
derstanding its nature. One such decay, A (3872) — > 
J/ipli IE EH established its charge parity to bo +1. 
In the molecular model, the radiative decays of the 
A (3872) occur through vector meson dominance (VMD) 
and light quark annihilation (LQA) [§[. The decay rate 
of A (3872) ~> JNry is dominated by VMD while for 
A (3872) —> ip'~f [l3j, it is mostly driven by LQA, im- 
plying that A" (3872) decay to 4>'l is highly suppressed 
compared to JNry |8|]. Recent results from the BaBar 
Collaboration Q show that #(A(3872) -> ip'j) is al- 
most three times that of £i(A(3872) — > J/ybj), which is 
inconsistent with a pure D*°D° molecular model, and 
can be interpreted as indicating a cc-D*°D° admixture 
[1, UK]. If the A (3872) is an admixture of x'd an d a 
molecular state, and its production and radiative decays 



are mainly due to its x'ci component, then the ip ,r y decay, 
a favored El transition of x'c\- should be significantly en- 
hanced compared to the J/ipj decay, which is "hindered" 
by poor wave function overlap [r|. 

In this Letter, we present new results on B — > 
(Xd, Xc2, A (3872)) A, where the X a, Xc2, A(3872) decays 
to J/V>7, and the A(3872) decays to V'7 [13 ■ These re- 
sults are obtained from the final data sample of 772 x 10 6 
BB events collected with the Belle detector [l8[ at the 
KEKB [lj| energy-asymmetric e + e~ collider operating at 
the T(45) resonance. The Belle detector is a large-solid- 
angle spectrometer which includes a silicon vertex de- 
tector, a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cherenkov counters (ACC), time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter (ECL) comprising CsI(Tl) crystals lo- 
cated inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. 

The J/ip meson is reconstructed in its decays to i + £~ 
[l = e or fx), and the ip' meson in its decays to i + £~ 
and J/ipir + TT~ . In the ip 1 — > e + e~ and J '/ip — > e + e~ de- 
cays, the four-momenta of all photons within 50 mrad of 
each of the original e + or e~ tracks are included in the in- 
variant mass calculation [hereafter denoted as M e + e -^], 
in order to reduce the radiative tail. The reconstructed 
invariant mass of the J ftp candidates is required to sat- 
isfy 2.95 GeV/c 2 < M e+e - (7) < 3.13 GeV/c 2 or 3.03 
GeV/c 2 < M M+Al - < 3.13 GeV/c 2 . In the ip' -> £+£~ 
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reconstruction, the invariant mass is restricted to the 
range 3.63 GeV/c 2 < M e+e - h) < 3.72 GcV/c 2 or 3.65 
GeV/c 2 < M 



< 3.72 GcV/c 2 . To reconstruct 
tp' — > J/ijm + Tr~ decays, AM = M#+£-„+„- — M t +i- 
should satisfy the condition 0.58 GeV/c 2 < AM < 0.60 
GeV/c 2 . In order to reduce the combinatorial back- 
ground due to low-momentum pions, the invariant mass 
of the two pions from the ip' decay, Jl^+j-, is required 
to be greater than 0.40 GeV/c 2 . A mass- and vertex- 
constrained fit is performed to all the selected J/ip and 
ip' candidates to improve their momentum resolution. 

The Xa,c2 and the A (3872) candidates are formed by 
combining the J j %p candidates with a photon. The pho- 
tons are reconstructed from energy depositions in the 
ECL and are required to have energies (in the lab frame) 
greater than 270 (470) MeV for Xd,c2 (A"(3872)) recon- 
struction. In a similar fashion, A (3872) candidates de- 
caying to f/ ; '7 are reconstructed by combining -0' candi- 
dates with 7 candidates with energies greater than 100 
MeV. 

Charged tracks are identified as pion or kaon candi- 
dates using information from the CDC (dE/dx), TOF, 
and ACC systems. The kaon identification efficiency 
is 88% while the probability of a pion misidentified as 
a kaon is 10%. The pions used in the reconstruction 
of the ip' in the J/ipTr + ir~ channel have an identifica- 
tion efficiency of 99% with a kaon to pion misidentifi- 
cation probability of 2%. Candidate Kg mesons are re- 
constructed by combining two oppositely charged tracks 
(with a pion mass assumed) with invariant mass lying be- 
tween [0.482,0.514] GeV/c 2 ; the selected candidates are 
required to satisfy the criteria given in detail in Ref . fiol ] . 

To reconstruct the B candidates, each J/V>7 or ip'^f 
system is combined with a kaon candidate. Two kine- 
matic variables are formed: the beam-constrained mass 



p* B 2 ) and the energy difference (AE 



E* B — E^ cam ). Here E^ cam is the run-dependent beam en- 
ergy, and E* B and p* B are the reconstructed energy and 
momentum, respectively, of the B meson candidates in 
the T(4S I ) center-of-mass (CM) frame. Candidates hav- 
ing Mbc > 5.27 GeV/c 2 and within a AE window of 
[-25,30] MeV for X d,c2 and [-30,35] MeV ([-20,20] 
MeV) for A(3872) -V J/V>7 (A(3872) -> V'7) are re- 
tained for further analysis. We extract the signal yield by 
performing an unbinncd extended maximum likelihood 
fit to the variable defined as M m-, — Mu + [2l[ , 



where is the world average mass 22j . In order to im- 
prove the resolution of M^, 7 , we scale the energy of the 
7 so that AE is equal to zero. 

To suppress continuum background, events having a 



ratio of the second to zeroth Fox- Wolfram moments [23 1 
i?2 > 0.5 are rejected. Large B — > tpX MC samples 
(corresponding to 50 times the data sample size used in 
this analysis) are used to study the background. To study 
the non-J/0 (non-0') background Mu sidebands in data, 
within [2.5-2.6] GeV/c 2 ([3.35-3.45] GeV/c 2 ) and [3.2-3.5] 



GcV/c 2 ([3.8-4.0] GcV/c 2 ), are used. 

For the (J/tpj)K channels, the background is primar- 
ily from B — > J/ipK* decays that do not peak in Mj/^. 
To reduce this background, we veto candidate photons 
from 7T° — > 77 by combining them with any other photon 
and then by rejecting both 7's in the pair if the ir° like- 
lihood is greater than 0.52. This likelihood is a function 
of the laboratory energy of the other photon, its polar 
angle and the invariant mass of the two-photon system, 
and is determined using MC study (24[. We also reject 
photon candidates with cos0i lc i > 0.76 (> 0.85) in the 
Xci,c2 (A(3872)) selection, where the helicity angle #hei 
is defined as the angle between the direction of the pho- 
ton and the direction opposite to the B momentum in the 
Xci,c2 (A(3872)) rest frame. Applying these criteria, the 
background is reduced by 86% (79%) with a signal loss of 
35% (30%) for the B X ciK (B A(3872)A) decay 
mode. For 1.3% of events with multiple candidates in 
B —> (J/ipj)K decay modes, we select the B candidate 



22] 



having M^ c closest to the nominal B mass 

A sum of two Gaussians is used to model the signal 
shapes of B — !> XciK and B — > Xc2K. The fraction of 
each Gaussian is fixed to the value obtained from MC 
simulated events. For B + — > XciK + the other shape pa- 
rameters are floated in the fit whereas for B + — > Xc2_K. + 
they are fixed using the mass difference (from Ref. [221 ] ) 
and the width difference (from MC simulations) between 
the Xci and Xc2- The non-peaking combinatorial back- 
ground component is modeled with a second-order poly- 
nomial. For the B° — > \ciKg and B° — > Xc2^s decay 
modes, the signal shape is fixed using the results from 
the charged B mode. 
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FIG. 1: M 
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Xcl,c2(- 



J/ip-y)K+ and (b) B u ->• Xd,c 2 (-> J/HWs decays. The 
curves show the signal (pink dot-dashed for \ c \ and red 
dashed for Xc2), and the background component (black dot- 
ted) as well as the overall fit (blue solid). The insets show a re- 
duced range of Afj/^, 7 and the contribution of the B — > x^K 
peak. 
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Figure 1 shows the fit to the Mj/^ distribution for 
B — > XciK and B — ¥ Xc2K decays in the range of 
[3.38, 3.70] GcV/c 2 . We observe the Xd in both B de- 
cay modes, and obtain 3.6 standard deviation (a) evi- 
dence for the Xc2 in the charged B decay mode. The 
statistical significance is defined as \J — 2 ln(£o / £ max ) 
where £ max (£o) denotes the likelihood value when the 
yield is allowed to vary (is set to zero). The system- 
atic uncertainty, which is described below, is included in 
the significance 25j. As no significant signal is found 
for B° — > Xc2K°, we determine a 90% confidence level 
(C.L.) upper limit (U.L.) on its branching fraction with 
a frequentist method that uses ensembles of pseudo- 
experiments. For a given signal yield, 1 0000 sets of signal 
and background events are generated according to their 
PDFs, and fits are performed. The U.L. is determined 
from the fraction of samples that give a yield larger than 
that of data. 

For the B X(3872)(-> J/ip"f)K decay mode, a 
sum of two Gaussians is also used to model the sig- 
nal PDF and the combinatorial background component 
is modeled by a first-order polynomial. To take into 
account small differences between the MC simulation 
and data, the signal PDF shapes are corrected for cal- 
ibration factors determined from the B + — > XciK + 
fit. Figure 2 shows the fit to the Mj/^ 7 distributions 
for B X(3872)K performed in the range [3.7, 4.1] 
GeV/c 2 . We find a clear signal for A(3872) -» J/ip-f 
in the charged decay B + — > X(3872)K + with a sig- 
nificance of 4.9(7 and measure the product branching 
fraction B(B+ — » X (3872) K+)B{X {3872) J/tp-y) = 
(1.78±£2(stat.) ± 0.12(syst.)) x 10~ 6 . We also give an 
U.L. on the branching fraction for the neutral B mode 
whose significance is 2.4c (Table HJ. We estimate the 
significance of the A(3872) — > J/ip'f signal by simultane- 
ously fitting the charged and the neutral B decay modes; 
we obtain a significance of 5.5a including systematics un- 
certainties. 
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FIG. 2: M 



j/i>~f 



distributions for (a) 



X(3872)(- 



J/ip-y)K+ and (b) B° -> X(3872)(-> J/ip^)K% decays. The 
curves show the signal (red dashed) and the background com- 
ponent (blue dotted) as well as the overall fit (blue solid). 



For the B — > (ip'^j)K decay mode, the background 



has a broad peaking structure, most of which is from 
B — > ip'K* decay mode. Here, since the 7's from 
A (3872) — > 1//7 have low energy (less than one third 
of the energy of the 7's coming from A (3872) — > J/ipj), 
the 7r°-veto and cos#hei selection result in more signal 
loss than background reduction. Instead, we combine 
the ip'K of the ip'jK candidates with any tt^ 1 or tt° 
candidate in the event. Three variables, namely AE' 

(= Efy+E* K . -E^J, M' hc (= ^E*l m -(p^+p* K ^) 

and the invariant mass of Kir {Mktt), arc used for this 
purpose. Events satisfying the criteria of 817 MeV/c 2 < 
M K tt < 967 MeV/c 2 , AE' within [-20,20] McV and 
M£ c > 5.27 GcV/c 2 , are identified as B -> ip'K* candi- 
dates and discarded. This results in the reduction of the 
background by 59% with a 22% loss of signal. For 15.4% 
of events with multiple candidates in B — > (ip'j)K decay 
modes, we select the B candidate having Mbc closest to 
the nominal B mass [22l |. 

The branching fraction for the B — > (ip'j)K mode 
is determined from a simultaneous fit performed to the 
two decay modes of the ip' . The background shape for 
B —> (ip'j)K has both a peaking and a non-peaking 
component. For the peaking component, the shape is 
estimated from a large sample of MC simulated events 
of ip'K and ip'K*, and their fractions are fixed using 
the branching fractions from Ref. [22[. The non-peaking 
background (combinatorial background) is parameter- 
ized by a threshold function (M^,^) 2 x cxp(a (M^,/ 7 — 
M Th ) + b (M^-Mrhf), where M Th = 3.725 GeV/c 2 . 
The ip mass data sidebands and large B — > ipX MC sam- 
ple (after removing B — > ip'K and B — > ip'K* decays) are 
used to estimate the parameters of the threshold func- 
tion. The shapes for both background components are 
fixed whereas their yields are allowed to float in the fit. 
The signal is described as a sum of two Gaussians and 
is fixed from MC study after applying calibration correc- 
tions (from B + — >• XciK + study) while its yield is allowed 
to vary in the fit. No significant bias is found in fitting 
ensembles of the simulated experiments containing the 
signal and background components. 

Figure 3 shows the results of the fit to the M^/ 7 dis- 
tribution for B -> X(3872)K. The fitted yields are 
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events) for B+ -> A(3872)A+ 
A(3872)Ag). Since there is no signifi- 



5.0 

(B 

cant signal in cither channel, we determine U.L.s of 
B(B+ -> A(3872)A+)Z3(A(3872) -> V" 7 ) (B(B° -+ 
A(3872)A°)£(A(3872) -> ip'~/)) as 3.45 x 10~ 6 (6.62 x 
10 -6 ) using the method described above. A completely 
independent analysis, with different selection criteria and 
a different fitting technique was performed on the same 
data sample; the results were found to be consistent with 
the results reported in this Letter. 

The branching fractions and the fit results are summa- 
rized in Table[l] Equal production of neutral and charged 
B meson pairs in the Y(4S) decay is assumed. Secondary 
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FIG. 3: M^/ 7 distributions for (a) B + -»■ A(3872) (-» 
ip'j)K + and (b) B° -» X(3872) (-»• ip'^K . The curves 
show the signal (red dashed for A (3872)) and the background 
component (pink dot-dashed for background from B — > ip'K* 
and B — > ip'K component, and black dotted for combinatorial 
background modeled by the threshold function) as well as the 
overall fit (blue solid). 



branching fractions used to calculate B are taken from 
Ref. fH. 



TABLE I: Corrected efficiency (e), signal yield (Y) from the 
fit, measured B or 90% C.L. upper limit (U.L.) for B — > 
Xa,c2K, B -> A(3872)(-> J/ipi)K and B -> A(3872)(^ 
ip'~/)K decay modes and significance (<S) with systematics 
included. B for B -> A (3872) A is the product B(B -» 
X(3872)A)S(A(3872) -> ^7). For B, the first (second) error 
is statistical (systematic). 



Decay e(%) Yield (F) 


Branching fraction 


S{o) 


B -) 


X ci(-» -7/^7)* 


B (xlO" 4 ) 




A+ 


14.8 2308^52 


4.94 ±0.11 ±0.33 


79 


K° 


13.2 542 ± 24 


3.7S±g;^ ± 0.33 


37 


B -s 


Xc2(^ J/ipl)K 


23 (xUT b ) 




A+ 


16.6 32.8i;o.2 


l.lll^ ± 0.09 


3.6 


A' 


14.4 2.8+3;g 


0.32i^±0.03 (< 1.5) 


0.7 


B - 


• A(3872)(^> J/ip-y)K 


B (xHT b ) 




A+ 


18.3 30.0t?.4 


1.781^ ±0.12 


4.9 


A 


14.5 5.7±|;| 


1.24+g;^ ±0.11 (< 2.4) 


2.4 


B - 


► X(3872)(^ i\)'i)K 


B (xHT 6 ) 




K+ 
K° 


14.7 5.o±ii-» 

10.8 1.5±|;| 


0.83ii;H ±0.44 (< 3 - 45 ) 
1.12±1;|5 ±0.57 (< 6.62) 


0.4 
0.3 



A correction for small differences in the signal detec- 
tion efficiency calculated from signal MC and data has 
been applied for the lepton (kaon/pion) identification 
requirement. Samples of J/ip — > i + (r and D* + —> 
D°(K~it + )it + decays are used to estimate the lepton 
identification correction and the kaon (pion) identifica- 
tion correction, respectively. The uncertainties on these 
corrections are included in the systematic error. The er- 
rors on the PDF shapes are obtained by varying all fixed 
parameters by ±ler and taking the change in the yield as 



the systematic error. To estimate the uncertainty arising 
from the fixed fractions of B — > ip'K and B — > ip'K* 
in the B — > (ip'"/)K background shape, we vary their 
branching fractions by ±lcr. The uncertainty due to the 
secondary branching fractions are similarly taken into ac- 
count. The uncertainty on the tracking efficiency and the 
number of recorded B meson pairs are estimated to be 
1.0% per track and 1.4%, respectively. The uncertainty 
on the photon identification is estimated to be 2.0% and 
3.0% for B (J/ipj)K and B {^'j)K, respectively. 
There is some possible efficiency difference of the selec- 
tions {E 1 , 7r°-veto and cos#hei) between data and MC. 
This difference in the B — » (Jfip'y)K study is estimated 
to be 3.0% using the B + — > XaK + sample. Due to the 
non-availability of a proper model to generate Xc2 m the 
Evtgen simulation [261 ] , and the ambiguity in the allowed 
X(3872) J PC values (1++ or 2 -+ ) [23], we generate Xc2 
and AT(3872) assuming them to be scalar, vector and ten- 
sor particles. We find that 4.0 % is the maximum possible 
difference in the efficiency and include it in the system- 
atic error. 

In summary, we observe A" (3872) — > J/ip"f in the 
B decays and present the most precise measurement 
to date of the product branching fraction B{B + — > 
X(3872)K+)B{X{3&72) -> J/ipy) = (1.78±g'j|± 0.12) x 
10~ 6 . We also report evidence for B — > Xc2K, and 
the ratio of B(B+ ->■ X c2K + )/B(B+ -» XdK+) is 
measured to be (2.25±£;^ ± 0.17)%. We find no ev- 
idence for A"(3872) — >• ip'j and give an U.L. on its 
branching fraction as well as the following limit i?(= 

< 2' 1 ( at 90% C - L ')- The X(3872) 
state may not have a large cc admixture with a D*°D° 
molecular component as was expected on the basis of the 
BaBar result 141. 
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